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ABSTRACT 


An interplanetary magnetic field, which should agree with observed cosmic ray and magnetic 
storm effects, has been suggested by Alfvén. This field is worked out more in detail in the present 
paper, and several cosmic ray orbits are calculated in this field. The results are briefly discussed 
and compared with cosmic ray observations. A beam field of 30 uw gauss at the earth’s orbit 
has been assumed here, and according to Venkatesan’s analysis of Forbush decreases the beam 
fields may vary from a few to about 100 mu gauss. 


Introduction 


Some time ago Alfvén (1956) suggested a type of interplanetary magnetic field, 
which might be compatible with observed cosmic ray and magnetic storm effects. 
It should be of great interest to calculate cosmic ray orbits in this field so that more 
detailed comparisons with observations in nature can be made. The numerical integra- 
tion has now been carried out at the Stockholm electronic computer BESK. 

We introduce a rectangular coordinate system with the origin in the centre of the 
sun, the z- and y-axis in the equatorial plane, where the beam is centered around the 
y-axis, and the z-axis towards the magnetic south pole (approximately the helio- 
graphic north pole). Spherical coordinates will also be used. The coordinate system is 


shown in Fig. 1. 


Beam direction 


Fig. 1. The coordinate system. $ 
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The magnetic field can be divided into different regions as shown in Fig. 2. It is 
assumed that there are no toroidal (g-) components of the field. The definitions of 
the regions are as follows. 


Region I. Polar regions at 0 <6 <6, and a — 6, <0 <a. The field is a dipole field. 

Region II. Beam region where an outward transport of flux takes place. By is 
assumed to be proportional to r-!. 6, <0<2—6,, —q@y<Y <p. 

Region III. Recovery region where an inward transport of magnetic flux takes 
place. The field recovers when the strongly perturbing beam has passed. By ~r-?. 
0, <0 <2 — 8s, Yo <Y <2 — Gp. 

Region IV. Transition regions between the beam region and the polar regions. 
6,<0<6, and m—60,<0<2—-6,, —% <9 <Q. 

Region V. Transition regions between the recovery region and the polar regions. 
6, <0 <6, and 7 —6,<0<2 —6,, 99 <p <2m —Q. 


Detailed calculations of the field 


For detailed calculations of the field we may use the fact that div B=0. We 
consider the sun as a homogeneously magnetized sphere, and the magnetic field 


inside the sphere may be By. The MKSA-unit system will be used. The sun’s radius 
is To. 
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. Regions I 

Polar regions with dipole field. The components of the field are 
B, = — B, cos 0 (rire), | 
Ba= — By sin 0 (r/ro)’, (1) 
5,0 


Region II 
Beam region. By is proportional to 7-1. Let us put 
Bo ae Kies. (2) 
Since B, vanishes according to our assumptions div B =0 gives 


é a 


r* = (° By) + (rsin 6) 56 (By sin 6) = 0. (3) 
This gives fBe=1, po er (By sin 6) dr. (4) 
aad sin 0 06 
"Oo 
B° = — B, cos 6 is the radial component of the field at the sun’s surface. Thus, in 


region II we get 


To \* ~—Y7, 
B,.= —'B, (‘) cos 0+ K, as cot 6, | 
Ba= = Ki As ' (5) 


B,=0. 


Region III 


Recovery region. Starting from the basic assumption that By ~r~*, we obtain in 
the same way as in region II 


2 
B,= — By (‘?) cos 0 + K,r~® cot 0 log — : | 


TO 
( 6 
Bo= — K,r”, | ( 


B,=0. 
Regions IV and V 


Transition regions. The field is according to our assumptions only determined by 


the following facts: 
(a) B, = 0. 
(b) div B=0. 
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Condition (b) makes it possible to connect the polar regions (I) with the equatorial 
regions (II and III) by means of field lines. For example, if we know that a certain 
field line enters the region at 7 = 12, 0 =9, at the border of an equatorial region, we 
can find the point r =r,, 9 =0, where this field line leaves the region. However, this 
is not enough to determine the field in detail. We need another assumption, and this 
may be made in such a way that it facilitates the calculations. We assume that the 
field lines are determined by 

sin 6 —sin 0, 
r ; 
sin J. — sin 0, 


(7) 


r=, + (1, 


Fig. 3 shows the field lines in regions IV and V qualitatively. Field lines through 
r=fro, 9=6, divide the transition regions into two different parts, a and }, 

(a) where the flux leaves the region through the cone  =6, (on the northern 
hemisphere), and 

(b) where the flux leaves the region through the solar surface. 

The field lines through r = ro, 6 =6, reach the cone 6 = 6, at r = 79, say, where 79 
is determined by the flux condition 


| Bo2ar sin 0,dr= — Byx712, (sin? 6, — sin? 6,) 
+O 
= flux through solar surface between 6, and 6,. (8) 


We get two different values of rj, one in region IV, 7p’, and one in region V, 7, 
corresponding to B, = — K, r~ (as in region II) and By = — K, r-? (as in region III), 
respectively. The equation (8) gives 


Byrd sin? 0, —sin? 6, 
2K, sin 0, ; 


IV 
(Ki =F) ae 


(9a) 


[Bord sin? 6, —sin? 6,| | 
Pad, snf, | 


r) =1o exp 


(9b) 


By calculating the flux through the cones 9 =6, between r =ro and r =71,, and 
6=06, between r =r, and r =r, we can get a connection between r, and r, i.e. between 
the points where a field line leaves and enters the transition regions [Va and Va 
(Fig. 3). At 6 =0, we take B, from (1). Thus, the flux through 6 = 6, is 


TY) 


= | By 2arsin d,dr= a Byr% sinea, (2-1), (10) 
: TO Fs 
fo) 

At 6 =, in the regions IT (IV) and III (V) we get the flux ¢, 

rv 

2 =| —K,r}2ar sin 0,dr= — 22K, (r3Y —719’) sin 6,, (11a) 
IV 
"O 
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Region I 


Dipole line 
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Fig. 3. Qualitative sketch of the magnetic 
field lines in the transition region. 


a 
py = —Kyr*2ar sinO,dr= — 22K, (log”y) sin 6, 
To 
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ry =179 exp 


{Bors sin? 0, (1 < “e)| 
| 2K, sin 0, 


Bor? in? 
=o exp | ae (sin 2 = a) re)! 


(ry,82) 


Region II or IL 


(11) 


(12a) 


(12b) 


The values in MKSA units of the different constants may in a representative case 


according to Alfvén (1956) be 
B, =1.2 x 10-3 Vs/m? = 12 gauss, 

' K, =0.65 x 10-6 ro Vs/m? ro (30 uw gauss at the earth’s orbit), 
K, =0.88 < 10-4 72, Vs/m? r3 (19 w gauss at the earth’s orbit), 


6, =60°, 
uP — 70°, 
PY = 15°. 


We define some new constants: 


oa 
_ Boro sin? 9, 


= —— = 136.7 ro 
“1 2K, sin 0, muro: 


(13a) 
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pe Ce 
Boro sin? 6, 


= = 5.442. (13b 
2K, sin6, 


X 


Introducing these constants in (12a) and (12b) we get 


Here 2 (1-2), (14) 
1 
ry =1) exp fe (1-29)! ; (14b) 
1 


Each field line in a certain meridional plane can be defined by the corresponding 1. 
If we want to calculate the magnetic field at the point (7,0,q) within the transition 
regions, we must first calculate the corresponding r,, i.e. we must know which field 
line goes through (r,0,~). This can be made by combining (7) and (14a) for region IV 
or (14b) for region V. We then obtain in region IV 


ae sin 9, — sin 6, p+ |e +4e4r6 


(sin 6, — sin 0) (sin 0 — sin 6,) \ (15) 
2 (sin 0,—sin 6) | 


(sin 6, — sin 6,)? 


sin # — sin 6, 


ak IV ‘BI U = SIN Uj) e 
where D==—r-+ (16. + &) cin Peetanedl | (15a) 
and in region V, rj is determined by 

a—bx=e i (16) 
eh _ sin #,—sin 6, r et, (16a) 


sin 9—sin 0, ro 


sin 0.—sin 8 ro 
neste : ae ee i6b 
sinf—sin@,ro 7 ’ ( ) 


ry 


a= (16 c) 


X27 


These equations will be used also in the regions IVb and Vb (Fig. 3) where the field — 
lines never reach the cone § = 6,. The values of 7, obtained there have no physical 
meaning but they are necessary for the calculations. 

After having computed r, we can calculate B in the following way. Consider a 
flux tube between the field lines r, and r, + dr,, and between two arbitrary meridional 


planes separated by an angle of one radian. The flux in this tube through the cones 
6 =6, and 6 =6 is determined by 


Bor sin 6 dr = By. 7, sin 0, dry. (17) 
Thus a sin 6, 7, dr, 
o= Bo, sinO r dr- ee 
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Fig. 4. The magnetic field lines in the beam. 


Bo, is determined by equation (1) with = 6,. In region IV we get from (7) and (14a) 


dr sin@,—sinf %?o sin#—sin 0, 


dri’ sin@,—sin6, (ri)? sin 0,—sin 6,’ ee) 
In region V (7) and (14b) give 
dr sn6@,—sin@ afro y sin #,—sin 6 
dry sin 0,—sin6,  (r¥)? ( ~ ? sin 0.— aa, ’ snd 
B, is determined by differentiating (7) since 
B,| By =dr/rd0 (20) 
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Fig. 5. The magnetic field lines outside the beam. 
fe—T; cos 8 
so that B,= Be (21) 


r sin 6,—sin 0,’ 
and this applies both to regions IV and V. 
The shape of the field lines is shown in Figs. 4 and 5. 
Calculation of orbits 


A modification of St6rmer’s method has been used for the orbit calculations. 
The velocity vector v can be divided in two parts 


Verve Vir (22) 


where v, is the component perpendicular to the magnetic field, and vy, is the component 
parallel to the field. 

When the integration has been carried up to the point r,, = (a,,Y,,;2,), Where the 
velocity is v,, then the next point r,,, is determined by 
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Payee tno Ata; (23) 


osin A « As 
tl oe ee 


Ar, =en(1—cosAa)+v Vall so (24) 


Oar 
with 


m 
0 =F (vx B) 


= radius of curvature of the motion perpendicular to the field. (25) 


As = step length of are along the orbit. 


Asv 
Ag=— 4, (26) 
Qo v 
m = relativistic mass of particle, 
e =electronic charge. 
Further we must compute V,,,,: 
sin A « 
Tae TESS mueti SS ae seacanalat (27) 


The modification of the Stormer method is that the magnetic field is not calculated 
at the point r, but at a point intermediate between r, and r,,,, except at the first 
step. This point is determined by 


Pn41/2 =Fn-4 “3 3 A Ini (28) 


This means that the step length can be rather long without loss of accuracy. The 
step Jength As must be chosen by considering both the gradient of the field and the 
angle Aw, which should not be too large. In the present orbit integrations 


As =1/32, (29) 


where 7 = distance from the sun. This means that B does not change by more than 
10 % in one step. If (29) would give Ax > 0.25 radians, then As is determined by (26) 
with Aw = 0.25. 


Selection of orbits 


The orbits calculated are proton orbits with energies 10°, 10!° and 10! eV. Each 
orbit is calculated backwards, starting on the earth (neglecting the earth’s field). 
The electric field in the beam has not been considered, because the effects of itare 
already known well enough, as will be explained later. The earth was situated at 
four different positions with respect to the beam (all of them in the equatorial plane) 
one in the centre of the beam (position A), one outside (opposite) the beam (position 
B) and two at the border of the beam (positions C and D). This is shown in Fig. 6. 
Since the beam rotates together with the sun in an anti-clockwise direction as seen 
in Fig. 6, position D corresponds to the time just before the beam overtakes the earth. 
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Fig. 6. The four positions of the earth relative to the beam. 


Therefore, this is referred to as “‘before” the beam. For the same reason position Cc 
is referred to as “‘after’’ the beam. The direction of impact on the earth was, still 
neglecting the geomagnetic field: 


v, = (0, 0, —v), 
Vo4=(+0/V2, 0, —v/V2), 
v3.5 = (0, +y/V2, —v/V2), 

Ve10 = (+2, 0, 0), 
Vga2— (0, £2, 0), 


V7,9,11,13 = (£/ /2, +v/V2, 0). 


Thus, 3 energies, 4 positions of the earth and 13 directions of impact on the earth 
give 156 orbits. Each orbit has an identification number and the explanation of these 
numbers is given in Table 1. 

In each group of 13 orbits (same energy and same position of the earth) the last 
8 are confined in the equatorial plane. The first one is perpendicular to this plane when 
impinging on the earth, and the 4 others make an angle of 45° with it. The directions 


Table 1. Explanation of the orbit numbers. 
pa ee ee ei ee ee 
Position of the earth (see Fig. 6) 
Particle energy - 
eV A B C D 
(in the beam) |(outside the beam)| (after the beam) |(before the beam) 


1 


10° 1-13 14-26 27-39 40-52 
101° 53-65 66-78 79-91 92-104 
101 105-117 118-130 131-143 144-156 


PO SE eS 4 
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Fig. 7. The identification numbers of the orbits refer to the 
directions of impact on the earth. No. 1 comes perpendicular 
to the equatorial plane, 2—5 are inclined 45° and 6-13 are 
in the equatorial plane. To each number should be added 13 
nm, where n = 0, 1, 2,..., 11. 


of the orbits with respect to the x- and y-axis (the beam is ejected in the y-direction) 
are indicated in Fig. 7. 
The probable errors in the calculations are discussed in the Appendix. 


Results 


The results are given in the Plates I-XLIII. Some of the orbits are not given 
there, since they do not contribute anything essential to our knowledge. 
The unit of length, LZ, in the plates is such that 


L=10.241r6. 


The reason for this is that 1024 = 21° is convenient to use in BESK. Hence, 1 A.U. = 
215 ro = 21.0 L. 

It is seen that the motion of the particles in trapped orbits is made up essentially 
of oscillations along the magnetic field lines from the northern to the southern hemis- 
phere and back again, plus a drift around the sun, clockwise if seen from above the 
northern hemisphere. There is also a radial drift at the edges of the beam because of 
the different magnetic fields inside and outside the beam. It should be noticed, that 
for r>136 ro = 13.3 L, B is smaller outside than inside the beam. For r< 13675 
the opposite is true. Therefore, it may happen that a particle which goes almost 
perpendicular through the equatorial plane is displaced from the sun when leaving 
the beam, as is the case for orbit 27. There the main radial drift is caused by a displace- 
ment at r<136 ro above the equatorial plane, where the field inside the beam is 
smaller than outside. 


Discussion 


From the results of the orbit calculations some conclusions may be drawn, which 
can be directly related to observations in nature. 

At first some general remarks should be made. The mathematical model which 
has been treated here is an idealisation and in nature the following deviations from 
this model affect the results. 


(a) There are electric fields in nature, in particular in the beam. The effects of 
this are, (i) a deceleration corresponding to the potential difference across the beam, 


17] 


L. BLOCK, Cosmic ray orbits 


and (ii) a radial drift outwards from the sun with the velocity v» = E/B (beam 
velocity). The deceleration of the particles produces Forbush decreases when the 
earth is in such a position that it can be hit by decelerated particles. The radial 
displacement of a particle passing the beam due to the electric field is equal to the 
product of v) and the time that it takes to pass the beam. This time is longer, the 
lower the energy is. In fact, it is inversely proportional to the kinetic energy per- 
pendicular to the field, as long as the radius of curvature is small compared to the 
width of the beam (W, < 10! eV, say). For higher energies the displacement 1s neglig- 
ible. For particles with W, = 10° eV at 1 A.U. from the sun the diplacement is about 
0.4 A.U. if vy, =2 x 108 m/s and if our assumptions about the magnetic field and the 
width of the beam are substantially correct. 

The effect (i) of the electric field (the deceleration of cosmic rays) is of decisive 
importance for intensity variations of cosmic rays. Hence, an important part of the 
application of the present orbits to observations in nature is to investigate whether 
cosmic rays impinging on the earth from a certain direction at a certain time have 
recently traversed electric field regions (mainly beam regions) or not. 


(b) The strength of the magnetic field may differ from that assumed here. In 
fact, Venkatesan (1957) has obtained values of the magnetic fields in the beams from 
studies of cosmic ray variations during magnetic storms. He has obtained values 
both from statistical means of many disturbances and from individual decreases. In 
the statistical means the intensity changes have been smoothed out and, hence, 
these result in much lower fields than the individual decreases. It seems reasonable to 
assume that the fields obtained from the individual decreases are more directly related 
to beam fields and, therefore, we will compare these with the assumed fields in the 
present model. For this we have to use Tables 6 and 7 in Venkatesan’s paper. These 
tables give a very wide range of values, varying from time to time. The magnetic 
fields obtained are inversely proportional to the assumed velocity of the beam. 
Venkatesan has given values for v, = 2 x 108, 108, 5 x 10° and 10° m/s. For2 x 106 m/s 
the fields obtained range from 1.45 x 10-* up to 78.4 x 10-® gauss as compared to 
our value of 30 x 10-* gauss. Now, 2 x 108 m/s is probably about the highest beam 
velocity, that ever occurs, and it is only reached at very strong storms. Hence, the 
lower field values, which refer to the weakest storms, should probably be connected 
with lower velocities, so that in reality the fields are not quite as low in these cases. 
It should also be noticed that Venkatesan has only obtained the magnetic field 
component perpendicular to the motion of the beam. The field may have other com- 
ponents, so the total field is larger. Nevertheless, it seems as if the beam fields 
assumed by Alfvén are rather high, so that they should be referred to rather heavy 
storms. 

In view of this we may say that the magnetic fields may differ from those assumed 
here by a factor 3, up or down, say, or in extreme cases by a factor 10, at least down. 
If the field is a third of what has been assumed, the orbits correspond to energies 


of roughly 1.6 x 108, 2.8 x 10° and 3.3 x 101° eV, and if it is 3 times larger the energies 
are about 4.5 x 109, 3 « 101° and 3 x 1011 eV. 


(c) The beams may not be symmetric around the equatorial plane, and they are cer- 
tainly not as regular, and the borders are not as well defined as assumed here. In par- 
ticular the surface currents connected with the sharp bends of the field lines at the bor- 
der of the transition regions do not exist in reality. However, the latter is likely to be less 
important for the orbits than the electric field and the strength of the magnetic field. 
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(d) The earth may not be situated in the equatorial plane but some distance from 
it, viz. up to 7°, which is the inclination of the ecliptic to the sun’s equator. The 
magnetic equator may also differ somewhat from the rotational equator. The beam 
may be ejected in such a direction that it does not hit the earth but passes over or 
under it. In such a case a Forbush decrease may occur without an associated magnetic 
storm. 


(e) Here it has been assumed that the beam is straight in the y-direction whereas 
in reality it is curved because of the sun’s rotation. It may very well be, that the 
field lines are curved in the same way too. We may say, that the magnetic meridional 
planes are curved. The direction of the curvature may change too, because if matter 
is moving inwards outside the beam this should result in a curvature in the opposite 
direction. 


(f) The exponent of the field has been assumed to be constant (=2) in the whole 
recovery region. It is likely that it varies, so that it is larger, the longer the time since 
a beam passed over the region in consideration. This means e.g. that the exponent 
should be larger in the recovery region near position D (Fig. 6) than near position (. 

(g) There may be no beam or several beams at the same time. 


These differences should be kept in mind when applying the results of the calcula- 
tions to observations in nature. 


1. The flare effect 


A remarkable flare effect occurred on Febr. 23, 1956. This has been studied among 
others by Meyer, Parker and Simpson (1956) and by Eckhartt (1957, 1958). The 
former assume a field free region between sun and earth so that the rays can easily 
move from the sun to the earth in the short time observed between the flare and the 
onset of the cosmic ray intensity increase. However, Eckhartt found that the direc- 
tion of impact of the rays on the earth was not from the sun, but from a beam which 
produced a magnetic storm two days later. The beam appeared to emit cosmic rays 
(see Fig. 8). It is difficult to explain this if the region between sun and earth is field 
free. However, it is also difficult to explain it in terms of the present field, since it 
implies an anti-clockwise drift of the rays around the sun. Contrary to this all orbits 
integrated so far on the whole obey a clockwise drift direction. It is true, that for 
many orbits the clockwise drift is frequently interrupted by anti-clockwise drifts 
(see e.g. orbits 82 and 96, Pl. XXX and XXX1I), but these are always so small that 
they cannot explain the effect observed by Eckhartt. Further, the anti-clockwise 
drift always takes place further from the equatorial plane than the earth is ever 
likely to be. This means that one would have to assume that the particles first 
drift anti-clockwise to a point where they can move along a field line towards the 
earth. However, then the rays impinging on the earth would appear to come from 
the sun. Hence, in this way the drift from the beam would be explained but not the 
observed direction of impact on the earth. 

Eckhartt (1957) has worked out an explanation in terms of scattering by the beam 
of cosmic rays from the sun. He shows how 15 GeV particles can move from the beam 
to the earth if the field outside the beam at the earth’s orbit is 1.9 x 10-® gauss. 
This is a rather weak field, only one tenth of the field assumed in the present paper 
(although not inconceivably weak, in view of Venkatesan’s result), The rays also 
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Fig. 8. The earth and the solar beam which is sup- 
posed to be responsible for the magnetic storm on 
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contain particles with lower energies, so actually one has to go down even more with 
the field. On the other hand, the field in the beam must be very strong in order to 
be able to scatter 15 GeV particles, at least as strong as assumed in this paper. This 
means a very strong gradient at the edge of the beam. 

There is an alternative possibility, viz. that the field lines are not confined in 
meridional planes but that they are curved, about as shown by the dashed line in 
Fig. 8. Then we do not have to assume very weak magnetic fields, because even very 
low energy particles can drift along the field lines, and the correct impact direction 
on the earth is automatically obtained for particles of all possible energies. 


2. Forbush decreases 


It has been pointed out long ago (Brunberg and Dattner, 1954) that most of the 
cosmic ray particles, which reach the earth just after a beam has passed by, must 
have traversed the beam where they were decelerated by the electric field as explained 
earlier. This would cause Forbush decreases. In agreement with this it has also 
been shown by Venkatesan (1957) that during high solar activity there is a strong 
correlation between K, and the negative time derivative of the cosmic ray intensity. 
Now it can be seen from the orbits 27-39 (Pl. X-XIV and XIX) and 79-91 (PI. 
XXVII-XXXI and XXXVI) that in general trapped orbits have to pass the beam. 
However, at 10" eV the orbits are not trapped and then not all of them have to pass 
the beam in order to reach the earth just after the beam has passed by. This can be 
seen from orbits No. 131-143 (Pl. XX XVII, XX XIX, XLII). The only orbits which 
definitely have to pass the beam are those reaching the earth from between the 1-hour 
and the 10-hour directions approximately. It is however, interesting to note, that all 
particles impinging from nearly the zenith direction at the magnetic poles of the earth 
must first traverse the beam. This can be seen from orbit 131, Pl. XX XVII, and it 
may imply that the Forbush decreases are more pronounced for higher energies 
at the poles than at lower latitudes, as seems to be the case according to measure- 
ments by Singer (1957) at Thule. 

The fact that the cutoff energy in the geomagnetic field is lower at the poles than 
at lower latitudes should also result in heavier Forbush decreases, since the lower 
energy particles suffer more intensity decreases because of the larger percentage 
of losses in the beam. This consequence of the electric field theory has also been pointed 
out by Singer (1957). In that connection he states, however, that a latitude depend- 
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ance of the Forbush decreases has not been observed. This seems to be in contradic- 
tion to his statement that the Forbush decreases are twice as large at the poles as 
in lower latitudes. If he means that the decreases do not depend on the latitude 
except near the poles, he may be right, but it seems as if sufficient data about this 
are not yet available. 

It is interesting to notice in Pl. XX XVII that even 10" eV cosmic rays impinging 
in the zenith direction at the poles never go very far from the solar equatorial plane. 
The maximum distance from the equatorial plane for these orbits is not more than 
0.5 A.U. except very far away (>2 A.U. from the sun) and the maximum solar 
latitude is about 25°. The same applies to many 10!! eV orbits inclined 45° to the 
equatorial plane (orbits 106-109, 119, 122, 133, 134, 146-148). 


3. Diurnal variation during Forbush decreases 


It has been remarked by Dattner and Venkatesan (1958), that the above mentioned 
fact, that a few orbits do not have to pass the beam in order to reach the earth after 
a storm, should cause a diurnal variation during Forbush decreases with a minimum 
around 6 h. They based their statement on orbits in the equatorial plane only. (Cf. or- 
bits 136-143, Pl. XLII). However, also the orbits No. 132-135, Pl. XX XIX (in- 
clined 45° with the equatorial plane) confirm their conclusion. 

Diurnal variations in general have previously been discussed by Brunberg and 
Dattner (1954, Sandstrém (1956a, b,) and by Dattner and Venkatesan (1958). The 
present model cannot directly account for most of these effects, since they depend 
on dynamic properties of the field and the cosmic ray gas, e.g. changes in the field 
with associated induced electric fields, diffusion due to density gradients and particle 
winds due to co-rotation with the sun. These effects may more successfully be studied 
by macroscopic magneto-hydrodynamic methods rather than by individual particle 
orbits. 


Conclusion 


The very brief discussion above indicates that the present model of the magnetic 
field may be a promising approach, which should be further investigated. 


APPENDIX 
Errors in the orbit calculations 


The errors in the orbit calculations are of two types: 


(i) Errors due to the limited capacity of BESK’s memory. 
(ii) Errors due to the numerical method of integration. These errors depend on 


the step length. 


The errors of type (i) are negligible in comparison to the other errors. The program 
for the machine is made such that the relative error in any number can never be 
larger than 10-7. By the relative error of a vector quantity we here mean the absolute 
error of this quantity divided by the modulus of the total vector. 

Thus we can focus our interest on the error due to the method of integration. 
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Refering to eq. (23), let us for a moment assume that Ar,, is parallel to the z;-axis, 
Then the error d(Aa;) is of interest only if 7+ 7. If there were errors only along the 
orbit, a plotted calculated orbit would never deviate from the true orbit. It is only 
the directional error 8 in Ar, which is important. Then 


it p i Se 
B= na V2 daa. he 


There will also be a directional error y in the new velocity vector v,,,, as calculated 
by (27). Omitting the index x +1 we find 


y=v? V2, {dv;}2, (31) 


so the total directional error per step will be 
6=Vp2+y2. (32) 


After n steps this will have increased to Vno. Thus, the total mean transverse error 
at the end of an orbit will be 


s 
A= l/ > {Vn As}?~NbAs=58, (33) 


where S=N As is the total length of the orbit. 

There are different cases with which we have to deal when determining 6. They 
depend on the five different regions of the field with which we are concerned and 
also whether the step length is determined by (29) or by Ax = 0.25. 

The calculations of 6 are lengthy but quite straight forward from (24)—-(27). In 
the two most important field regions the following approximate results can be given. 

Beam region: 
64, ~10-*(W +1) if Aa=0.25 
612 ¥ 1077/(W +1) if As=r/32 


Recovery region: 


3p, ¥10-°— (W+1) if Ax=0.25 


TO 
r 5 : 
Oo9 epety.4 al if As=r/32 


Here W is the energy in GeV. W +1 is a factor which is proportional to the mass 
of the protons, which happens to have a rest mass close to 1 GeV. At 1 A.U. from 
the sun the limit, where Ax =0.25 gives the same step length as As = 17/32, is at 
about W =10 GeV. Hence, as an average Ax =0.25 is used for W <10 GeV and 
As =1r/32 for W>10 GeV. This means that the maximum errors occur for about 
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10 GeV. Thus for r = 200 r5 ~1 A.U. and for 10 GeV 6 is of the order of 10-3 or 
2 x 10%. Since the orbits are about 1000 ro long, the errors will be one or a few solar 
radii according to (33). 

In the case of orbits entirely in the equatorial plane the true orbits are well known, 
so a comparison with the computed orbits is easy to make. The results are typed 
by BESK with 4 digits, where one unit in the last digit corresponds to 0.1 ro. Thus, 
the smallest deviation from the true orbits, which can be detected, is 0.1 ro. Of 32 
orbits calculated in the equatorial plane at 10° eV only one deviated as much as this, 
and in this case the deviation was 0.5 ro towards the sun and no deviation in the 
z-direction. This was for orbit No. 6. Here the applicable angular deviation is 5,, ~ 
2 x 104 so that the mean deviation after 10? ro is of the order of 0.2 ro. 

The largest errors should occur for 101° eV. Of the 32 orbits in the equatorial plane 
for this energy two showed no detectable deviation, namely No. 84 and 91. The 
mean error for orbits at 10! eV should be of the order 1-2.5 ro. Thus, the orbits 84 
and 91 deviate by less than one tenth of the probable deviation. No orbits show 
any detectable deviation when half of the integration is performed. After this they 
deviate more, sometimes even more than 10 76 but on the whole our estimate of 
the errors seems to be confirmed. 

For 101! eV no deviations have been detected. Thus, it is seen that for 101° eV 
the errors are so large that 10° 7, is about the longest orbit that can be safely integrated 
with this method. 


The author wishes to express his sincere thanks to Professor H. Alfvén, who has proposed 
this work, and to Dr. N. Herlofson whose great knowledge about orbits in magnetic fields has 
been of very good help in many discussions. The investigation has been made possible by grants 
from Statens Naturvetenskapliga Forskningsrad. 
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Plates I-XLIII 


The position of the earth is indicated by an arrow, except in a few cases where this would only 
confuse. The start and end points of each orbit are indicated by rings. Except in Plates XLII and 
XLIII (see below) dashed lines indicate approximate or simplified orbits, e.g. drift orbits. Examples 
of this are 10® eV orbits which oscillate along the magnetic field lines at the same time as they 
drift in the —q-direction. In this case the dashed lines indicate the approximate amplitude of the 
oscillations. 

Many 101° eV orbits are so complicated that it would be impossible to trace the orbits if they 
were plotted in their whole length. In these cases a part in the middle has been omitted. This 
omitted part may contain a few oscillations along the magnetic field lines. For orbit 19 only one 
point (indicated by a ring) on each loop of the trochoid has been given. The distance between two 
consecutive rings indicates the drift per revolution. 

The dashed lines in Plates XLII and XLIII indicate orbits not directly integrated but copied 
from corresponding orbits with the earth in position B. 
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